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We study the finite-size scaling of the roughness of signals in systems displaying Gaussian 1/f 
power spectra. It is found that one of the extreme value distributions (Gumbel distribution) emerges 
as the scaling function when the boundary conditions are periodic. We provide a realistic example 
of periodic 1/f noise, and demonstrate by simulations that the Gumbel distribution is a good 
approximation for the case of nonperiodic boundary conditions as well. Experiments on voltage 
fluctuations in GaAs films are analyzed and excellent agreement is found with the theory. 
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It was about seventy years ago that a current-carrying 
resistor was first observed to exhibit voltage fluctuations 
with a power spectra nearly proportional to the inverse of 
the frequency . Since then it has been shown that such 
an 1// noise is present in an extraordinary variety of phe- 
nomena J3|, examples being the light emission of white- 
dwarfs [g||7 the flow of sand through hourglass [Q, the 
ionic current fluctuations in membranes Q , and the num- 
ber of stocks traded daily Q. Most of 1/f fluctuations 
encountered are Gaussian, although non-Gaussian cases 
are known . The universality of this scale invariant 
phenomenon led to suggestions that part of the explana- 
tions should come from a generic underlying mechanism. 
Despite a large body of works, however, such a mecha- 
nism has not been discovered yet. 

In this letter, we establish a connection between 1/f 
noise and extreme statistics that may provide a new an- 
gle at the generic aspect of the phenomena. Namely, we 
shall show that Gaussian 1/f power spectra in periodic 
systems imply that the distribution of the fluctuations in 
the finite-size "width" of the signal is one of the extreme 
value distributions, the Gumbel distribution [Q. 

As we shall see, the Gumbel distribution emerges as a 
finite-size scaling function in the above connection. Thus 
the result can also be viewed as an interesting contri- 
bution to the gallery of nonequilibrium scaling functions 
that can be effectively used in investigating far from equi- 
librium processes. Indeed, imagine that a distribution 
function is measured in experiments (or simulations). 
Comparing this function with those in the theoretically 
built gallery (note the absence of fitting parameters) one 
can identify the relevant features of the underlying dy- 
namics in the experimental (or model) system. Such an 
approach was initiated in connection with surface growth 
problems |ll| , and the results have been used to estab- 
lish universality classes in rather diverse processes such 
as massively parallel algorithms [|l2| and the interface 
dynamics in the d — 2 Fisher equation ||l^. This line 
of reasoning ||l^ has also led to a connection between 
the dissipation fluctuations in a turbulence experiment 
Ilia and the interface fluctuations in the d = 2 Edwards- 



Wilkinson model [XY model) ||T^-[T^ and, furthermore, 
it helped in a parameter-free analysis of the upper critical 
dimension of the Kardar-Parisi-Zhang equation p9| . 

The derivation of the Gumbel distribution follows the 
steps of a similar calculation for the width distribution of 
random- walk interfaces ^ . Let the time evolution of the 
physical quantity of interest in the interval < t < T 
be given by h(t). This time series is equivalent to a 
surface configuration with h{t) being the height of the 
surface over a d = 1 dimensional substrate of length T 
with t being the coordinate along the substrate. The 
quantity of interest is the mean-square fluctuations of 
the surface (also called roughness or width-square of the 
surface) given by 



W2{h) = [h{t)-h] , 

where over-bar denotes average of a function over t, 



(1) 



(2) 



Let us assume now that the path probability of a given 
time series h{t) is known V[h{t)] ^ exp[— Then 
the probability distribution of the surface fluctuations, 
P{'W2), can be expressed as a path integral 

P{w2) = j Vh{t) 6 {w2 ~[h^ ]) V[h{t)] . (3) 

We shall restrict the above functional integral to periodic 
paths h{t) = h{t + T) and, in order to keep V normaliz- 
able, the integration is carried out with h kept fixed. 

The next step is to introduce the generating function 
for the moments of P{w2)'- 



G{s) 



dyP{y)e''y . 



(4) 



Substituting 'P[h{t)] ^ exp[— iS'[/i]] into (||) and evalu- 
ating the integral (|^), we find the following functional 
integral 
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G{s) =U Vh{t) cxp -S[h] -s{h^ - h ) 



(5) 



where A/" is a normalization constant to ensure G(0) = 1. 

The key question now is how to choose S in the prob- 
abihty density functional for 1// noise. The mathemat- 
ical representation of 1// noise has been pioneered by 
Mandelbrot and Van Ness , and since then a quite in- 
tricate theory (involving e.g. fractional derivatives) has 
emerged Due to the periodicity imposed on h(t), 

however, a simple form for S can be given in our case. 
Namely, we shall consider a "perfect" Gaussian 1// noise 
with the spectrum being linear for all frequencies, i.e. the 
following action is assumed 



S = a ^ |n||c„p = 2cr^n|c„| 



(6) 



n=l 



where cr is a parameter setting the effective surface ten- 
sion and the c„-s are the Fourier coefficients of the signal 



{N-l)/2 

h{t) = ^ c„e 

ri=-(Af-l)/2 



2iTint/T 



(7) 



Here h(t) is given on iV equidistant points {t ~ kAt, T = 
A^Ai), and we introduced the notation L = (N — l)/2 
with TV assumed to be odd. 

Using the functional integral can be written as 



G{s) = AT / V[c] exp 



E2( 



crn -t- s c„ 



(8) 



The integrals over the real and imaginary parts of c„ 
(n = 1,2,...L) yield a simple form for the generating 
function once we have used the condition G(0) — 1 to 
determine the normalization constant A/" in m): 



G{s)=Y[(l + — 



(9) 



The moments of P(w2) can now be calculated and, in 
particular, one finds that the average of W2 diverges for 
large L as 



[W2) = - 



dG 
ds 



1 ^ 



1 



InL 



(10) 



with 7 — 0.577... being the Euler constant. On the other 
hand, the fluctuations of W2 are finite 



{W2) = ^ 



TT 

7i 



(11) 



This means that the scaling variable 1^2/(^2) usually con- 
sidered [|| in cases when -^/(i 



{W2) is not the 

best choice since the distribution function reduces to a 



{W2 



delta function. Just as in case of the fluctuations of the 
d — 2 EW interface jl^ , the nontrivial shape underlying 
the delta function can be made visible by introducing 



W2 



{W2) 



{W2) 



(12) 



Then the L ^ 00 limit of the inverse Laplace transform 
of G(s) yields P{w2) in the following scaling form 



^x)^^{w',)-{w2)^P{w2)^ I l^re-nTri- 

— 100 

(13) 

The infinite product in (^ is equal to e('''/")r(l -I- s/a), 
so the inverse Laplace transform can be evaluated us- 
ing Euler's integral formula for the T function, and one 
obtains 



$(a:) = ae 



-(ax+'y) — e 



-(0X+7) 



(14) 



This scaling function, shown in Fig. 0, is one of the cen- 
tral results of our paper. In ^{x) one recognizes the 
Gumbel distribution [Q which is one of the three limit- 
ing forms of extreme value statistics. 

Extreme statistics has been studied in many contexts 
and the Gumbel distribution emerges frequently ||2^. A 
recent example in connection with surfaces is the study of 
the scaling behavior of the growth of the maximal relative 
height of a surface |24|. Since in most of the these stud- 
ies an extreme property is investigated, it is not entirely 
surprising to see the Gumbel distribution appearing. For 
1/f noise, however, this is not the case. A simple quan- 
tity such as the mean-square fiuctuations (roughness of 
the interface) is distributed according to extreme (Gum- 
bel) statistics. Although we do not see a simple physi- 
cal reason that necessitates this mathematical result, we 
speculate that it may be a key feature that underlies a 
unified treatment of systems displaying 1/f noise. 

When trying to compare the scaling function $(a;) with 
experimental results the following problem arises. An 
experimental signal is analized by moving a window of 
length L and building a histogram from the values of W2 
computed for the windows. The problem now is that 
the boundary conditions for the windows are not peri- 
odic and it is known that the boundary conditions affect 
the scaling functions . Thus two questions should be 
answered. First, is there a physical system with an effec- 
tive action (^ of the 1/f noise where periodic boundary 
conditions are realized? Second, how sensitive is $(a;) to 
the boundary conditions. 

Let us begin with the first question by showing an ex- 
ample of such a system. We consider the steady state 
fiuctuations of a d = 2 Edwards- Wilkinson (EW) sur- 
face with the substrate taken to be an infinite plane. We 
draw a circle of radius R on the substrate and compute 
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the probability density functional of a height configura- 
tion ho{ip) over this circle (parametrized hy < ip < 27r). 
We shall find that the action in this functional is equal 
to that of the Gaussian 1// noise given by (^). 
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FIG. 1. The analytical form ( p^ ) of the roughness distri- 
butions for a 1// signal in the periodic case (A'^ = oo, solid 
line) as compared to the numerical result for the nonperiodic 
case for several values of A'^ varying from 64 to 32768 (dashed 
lines). Inset shows the same curves on semilog scale. 

In the d = 2 EW model, the steady state fluctuations 
of the height of a surface h{x) are characterized by the 
following free-field action 

S2[h{x)] = ly \^h{x)\^d'x, (15) 

where the integral extends over a plane and the bound- 
ary conditions at infinity are fixed at h(oo) — 0. Working 
with polar coordinates (r, if) and specifying the height 
function on the r = R circle by h(r = R,ip) = ho{ip), the 
probability distribution functional of ho{ip) becomes 

P[ho{ip)] = I Vh[x) e-^^I''^^)! . (16) 

J h{S)\R=ho(ip} 

Here the functional integration is carried out with the 
boundary condition fixed at r = i? while Z is the same 
integral without restriction on the circle (thus P[/io(v)] 
is normalized to 1). Since the action is purely quadratic 
in the field, we have 

P[hoiip)] = Z-ie^^^[''<^(^)l / Vh{x) e-5^[''(=^)l, (17) 

Jh\B=0 

where hc{x) is the "classical" solution that extremizes 
the action with fixed boundary conditions, and the inte- 
gration variable h{x) = h{x) — hf.{x) . Since the integral 
does not depend on the boundary values, we can absorb 
it into the normalization coefficient, so the term of inter- 
est is the extremal (classical) action. The classical field 
satisfies the Laplace equation with fixed boundary con- 
dition, that is we have two Dirichlet problems, for both 



inside {h^ ) and outside (/i+) of the circle with a common 
boundary condition at r = i?: 

A/j±(f) = 0, hfU^hoiifi). (18) 

Let the solutions of the above problems be hf (r, ip) . Sub- 
stituting them into the action ( [l5| ) and applying Gauss' 
theorem, we can express the classical action S2[hc\ as 
functional of ho[Lp). Denoting this functional by S'c[/io], 
we have 

Sc[ho\ ^-J- J ^V^'^oif) [drh-{r,(p)\j^- drh+{r,(p)\^] . 

(19) 

The solutions of the Dirichlet problems are given by 

/.±(r,^) = -= J2 [-) hn,ne^-^, (20) 

where hn^n are the Fourier coefficients of ho(ip) satisfy- 
ing hfj^n = h*n _^. Substituting (|2^) into (|l^ yields the 
classical action for the height fiuctuations on the circle 

oo 

5c[/io(<^)] =2(T^n|/iK,„|2. (21) 

n=l 

Comparing ( ^ and (^, one finds that Sc[ho{'p)] is the 
action of a perfect Gaussian 1// noise. Since ho{ip) is a 
periodic function, we have thus indeed obtained a physi- 
cal realization of a periodic signal with 1// noise. 

Let us now turn to the question of how sensitive <i>(x) 
to changes in the boundary conditions. Analytically this 
turns out to be a hard problem and only numerical calcu- 
lations were performed. First, a periodic series of length 
N having 1 // power spectrum was produced using an ap- 
propriately filtered Gaussian white noise. Next, the sig- 
nal was divided into non-overlapping segments of length 
N and having determined the W2-s, the histogram of 
W2 and then the scaling function was built (using the 
same normalization as for the periodic boundary condi- 
tion case). In order to obtain satisfactory precision we 
used N — 2^** and averaged over 200 realizations of the 
periodic signals. 

For N — N the result for the periodic case is recov- 
ered while, in the N <^ N limit, one finds that $(x) is 
independent of the size of the segments in a wide range 
of N values. The ^{x) obtained in this N <^ N limit 
will be considered as the scaling function for nonperiodic 
(or "experimental") boundary conditions. As one can 
see from Fig. |l|, the distribution of nonperiodic signals 
deviates from that of the periodic case (Gumbel). The 
deviations, however, are small and mainly concentrated 
around the maximum of the function. 

We now consider the case of voltage fiuctuations in 
semiconductor films. The experiment was made by A.V. 
Yakimov and F.N. Hooge ||2^]. They considered n-type 



3 



epitaxial GaAs films grown by molecular beam epitaxy. 
A noise-free current passed longitudinally through the 
film and other contacts were used as voltage probes both 
for longitudinal and transverse directions. Several time 
series with typically 163840 points were obtained and 
the power spectrum was found to exhibit 1// behavior 
roughly over two decades of frequencies. We have re- 
analyzed these data by dividing the signal into segments 
of length N — 32, 64, 128, 256, and computing the dis- 
tributions of the roughness for different N. The results 
are displayed on Fig. ^. One can see that the experimen- 
tal data fit well with the theoretical curves (note that no 
fits are used in collapsing these functions) and that their 
precision is not good enough to distinguish between the 
periodic and nonperiodic cases. 




FIG. 2. Roughness distribution in the experiments on volt- 
age fluctuations in resistors calculated for A'' = 32 64 (x), 
128 (*), and 256 (□), compared to the analytical and numer- 
ical results shown in Fig. 0. Inset shows the same curves on 
semilog scale. 

As a final remark, let us note that the apparent ubiq- 
uity of the 1// noise is partly the result of loose terminol- 
ogy. Systems whose power spectrum is of the form 1//" 
with a close to 1 are also said to exhibit 1// noise. The 
approach we used for the treatment of the pure 1// noise 
can be extended to these systems and the dependence of 
the scaling function on a can be determined. A 

detailed discussion of this more general problem will be 
given in a separate paper |2^ . 
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